Heavy ion collisions in the used nucleon model 
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Abstract 

It is shown that recently proposed by R.J. Glauber the used nucleon model 
combined with the assumption that the nucleon consists of two constituents (a 



quark and a diquark) describes well the PHOBOS data on particle production 
at midrapidity. 
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1. In this note we investigate the consequences of the used nucleon model [1] in 
the context of the PHOBOS data on particle production at midrapidity [2-4]. This 
model was formulated as a simple generalization of the wounded nucleon model [5] 
which fails to describe correctly the observed particle multiplicities. 

Contrary to the wounded nucleon model, we assume that the number of parti- 
cles produced from one wounded nucleon does depend on the number of inelastic 
collisions this nucleon underwent. To formulate this model in more detail let us, for 
a while, consider nucleon-nucleus (mass number A) collision. For the incident nu- 
cleon contribution we have: the first inelastic collision produces n particles, where 
2n = (dn/dr))^^ is the average multiplicity at midrapidity in a single proton- 
proton collision. Following [1] we assume that the second collision produces a frac- 
tion /i of n, the third one /i 2 , and so on. The number of particles produced (from 
the incident nucleon) after k collisions is: 

1 ^ 

n k = n(l + /i + ...+ n k ~ l ) =n—^-. (1) 

1 — /i 
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From (1) it is seen that the used nucleoli model gives a natural way of interpo- 
lating between the two limits: for fi = we arrive at the wounded nucleon model [5] 
i.e. there is no difference whether a nucleon is hit once or several times, n& = n; for 
— 1 the incident nucleon contribution is proportional to the number of collisions, 

rik = nk. 

The probability Pk that the incident nucleon underwent exactly k inelastic col- 
lisions is given by a standard formula: 

P k (b) = [a NN T NN (b)] k [1 - a NN T NN (b)] A - k , (2) 

where T NN is given by: 

T NN (b) = — [ a NN (s) T(b - s)d 2 s. (3) 

&NN J 

Here ctjvjv and <Jnn(s) are the total and differential inelastic proton-proton cross- 
sections, respectively. T(s) is the thickness function (normalized to unity). 

The average number of particles produced from the incident nucleon is given by: 
E fc PkUk = 73^ { 1 " [1 ~ (1 ~ ^ a *NT NN (b)] A } ■ (4) 

Consequently, the average multiplicity at midrapidity iV = (dN / dr])^^ for the 
symmetric nucleus-nucleus collision at a given impact parameter b one obtain 

N(b) = 2A .^ U [T(b - s) {l - [1 - (1 - /x) a NN T NN (s)] A \ d 2 s, (5) 
a AA (b) 1 - /i J I J 

where aAA(p) is the inelastic differential nucleus- nucleus cross-section. 1 

2. The PHOBOS data are presented versus the number of wounded nucleons in 
both colliding nuclei W, given by [5] 

2A C 

W{ P) = m / T ( b ~ s ) i 1 ~ t 1 " °nnT nn (s)] a } d 2 s, (6) 

vaa{o) J 

Here and in (5) for the nuclear density we use the standard Woods-Saxon formula 
with the nuclear radius Ra u = 6.37 fm and d = 0.54 fm. We assume the differential 
proton-proton cross-section aNN (s) to be in a Gaussian form: 2 



a NN {s)= 1 e- s2 l"\ (7) 



1 In case of Au-Au collisions, using the optical approximation, we have verified that <JAA{b) = 1, 
for b < 14 fm. 

2 We performed the full calculation also in case of point-like interaction approximation, 
&nn{s) = <tnn5 2 (s)- We have observed that centrality dependence is slightly worse, indicating 
that the more complete treatment is indeed needed. 
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where x 2 = <Jnn/{^i) and 7 = 0.92 [6, 7]. 

We have observed that the used nucleoli model gives a good description of the 
RHIC Au-Au data with n ~ 0.41. 

The comparison of the used nucleon model (/x = 0.41) with the PHOBOS data [2- 
4] on the average multiplicity at midrapidity per one wounded nucleon as a function 
of W is shown in Fig. 1. Also shown are inelastic pp data (points at W = 2), as 
quoted in [3,4]. 

The inelastic proton-proton cross- sect ions needed for this calculation were taken 
as a NN = 32 mb, 36 mb, 41 mb and 42 mb at the RHIC energies y/s = 19.6, 62.4, 
130 and 200 GeV, respectively. 
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Figure 1: The predictions of the used nucleon model (for b < 14, corresponding to 
W > 5) with fx = 0.41 compared with the PHOBOS data [2-4]. The shaded areas 
reflect the inaccuracy in the inelastic pp data. 

3. To interpret the parameter \i let us consider two simple production scenarios. 

(i) Wounded quark model [8]. It is assumed that the nucleon consists of three 
constituent quarks and particle production from these quarks is independent on the 
number of interactions they underwent. Furthermore, we assume that the number 
of produced particles after k collisions is proportional to the number of wounded 
quarks, that is: 

(9) 

n w(i) ' { } 

where is the average number of wounded quarks per nucleon in a single inelas- 
tic proton-proton collision, is the average number of wounded quarks in the 
nucleon which underwent k inelastic collisions. The latter may be calculated by a 
straightforward counting of probabilities: 

w ( k q) = 3[ Pq + Pq (l - Pq ) + ... + Pq (l - Pgf- 1 } ■ (9) 
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Here p q is the probability for a quark to interact in a single proton-proton collision 
(note that p q = w^/3). Taking (8) and (9) into account we obtain: 

n k = n [1 + (1 - p q ) + ... + (1 - p q ) k - 1 ] = n l ~^~ Pq)k . (10) 

Pq 

Comparing (10) with (1) we find that fj, — 1 — p q . 

It turns out that at the RHIC energies approximately 1.2 [9], giving 

pi 0.6. We checked that this value leads to significantly larger multiplicities than 
actually observed [10]. 3 

(ii) Wounded quark-diquark model [7] . We assume that the nucleon is composed 
of two constituents (a quark and a diquark). As above we assume that particle 
production from these constituents is independent on the number of interactions they 
underwent, and both constituents produce the same number of particles. Performing 
analogous calculations as above we obtain: 

l-(l-p g ) k , 1-(1- Pd ) k nu 
n k = n ■ + n ■ , (11) 

Pq +Pd Pq+ Pd 

were p q and pd are the probabilities for a quark and a diquark to interact in a single 
proton-proton collision, respectively. 

As it is not clear what the diquark is, however, we may consider two different 
possibilities: 

(a) p d = p q , it corresponds to the situation where both constituents are similar. 4 
In this case (11) reduces to (1) with /i — 1 — p q . Assuming the average number 
of wounded constituents in a single proton-proton collision w^ q+d ^ = p q + Pd to be 
1.18 ~ 1.19 [7] we obtain the proper value /i ^ 0.41. 

(b) pd = 2p q , where diquark is rather large, comparable to the size of the proton 
(this relation was obtained in [7]). Now it is not possible to use formula (5), however, 
performing analogous calculations leading to (5) we obtain: 

N(b) = aUb) J**> J T (b-s){l-[l- Pq a NN T NN (s)] A + (12) 
+ 1 - [1 - Pd<?NNT NN (s)] A | d 2 s. 

We have checked that this formula withp^ = 2p q and = 1.185 [7] gives the 

results which practically do not differ (less than 5 %) than those presented in Fig. 
1. Thus, we may conclude that the model is almost independent on the details of 
the diquark 5 , and the only thing that really matters is the number of constituents. 

3 In order to obtain the value of fi ~ 0.41 we would have to assume that rj 1.75. This 
number, however is difficult to justify. 

4 It may mean that both constituents are of the same size. 

5 We also checked for other choices between pd = p q and pd = 2p q . 
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4. In conclusion, we have shown that the used nucleon model combined with the 
assumption that the nucleon consists of two constituents (a quark and a diquark) 
naturally describes the PHOBOS data on particle production at midrapidity. 
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